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ABSTRACT 

Context. Kepler produces a large amount of data used for asteroseismological analyses, particularly of solar-like stars and red giants. 
The mode amplitudes observed in the Kepler spectral band have to be converted into bolometric amplitudes to be compared to 
models. 

Aims. We give a simple bolometric correction for the amplitudes of radial modes observed with Kepler, as well as the relative 
visibilities of non-radial modes. 

Methods. We numerically compute the bolometric correction CK-boi and mode visibilities for different effective temperatures T^g 
within the range 4000-7500 K, using a similar approach to a recent one from the literature. 

Results. We derive a law for the correction to bolometric values: CK-boi = 1 + ai(7^ciT - T„) + a2(7^ctf - T„)^, with = 5934 K, = 
1.349 X 10""* K"' , and 02 = -3. 120 x 10"' K"- or, alternatively, as the power law CK-boi = (T^ff/Tof with a = 0.80. We give tabulated 
values for the mode visibilities based on limb-darkening (LD), computed from ATLAS9 model atmospheres for T^-ff e [4000, 7500] K, 
log^ 6 [2.5, 4.5], and [M/H] e [-1.0, +1.0]. We show that using LD profiles already integrated over the spectral band provides quick 
and good approximations for visibilities. We point out the limits of these classical visibility estimations. 

Key words. Asteroseismology - Stars: atmospheres - Stars: solar-type 
1. Introduction ening (LD), which depends on the observed wavelengths and 
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„ , , . ^ . , . . , r77 — ; n u^.^ . requests a knowledge of stellar atmospheres. 

Whde the NASA Kepler mission (e.g. jKocLetai; » is this note, we first propose in Lt. |2] a simple bolomet- 

. dedicated to exoplanet finding, the high-quality photometric correction for radial modes observed with Kepler, following 

^ data provided by the instrument are perfecdy suited for aster- approach similar to the one developed for the Conv ection, 

^ oseismology, specially for the study of s olar-hke oscil l ations Rotation, a nd planetary Transits (CoRoT, Ba dinetal.1 12OO6I1 

in main-seqiie ngg and red giant stars (e g. jCha plin et al.||2010t ^^^^5^^^ ^ iMichel etalJ (12009. hereafter M09). In Sect. |3] we 

| Beddingetal| 2010). For solar-like oscillations, the chai-acteris- visibihties of non-radial modes for Kepler, using LD 

. „ I'? °f '^'l^ frequencies, amplitudes, or ^.^^j^^ computed from ATLAS9 model atmospheres in a modi- 

^ . lifetimes are d etermned by fittmg the oscillation power spec- ^^^.^.^^ convection treatment, and we discuss the use 

^ trum (see, e.g.,| Ballot|20ia and references therein for a descrip- ^D laws integrated over spectral bands. Finally, we discuss in 

tion of these techniques). The apparent amphtudes of modes re- ^^^^ ^^^^ ^^^^^ correcting factors, and point out some of 

" ■ ■ covered in the power spectrum - or the light curve - depend on limits 
the spectral response of the instrument. It is therefore fundamen- 
tal to recover the bolometric amplitudes of modes to be able, for 

example, to compare the observed amphtudes to theoretical pre- 2. Bolometric Correction for radial modeS 

dictions. 

Moreover, owing to cancellation eff^ects, the apparent ampli- ^e derive in this part the bolometric coiTection for radial modes, 

tude of a mode depends on its degree {. Modes with degrees This correction is the factor to apply to the amplitude of a radial 

^ > 3 are generally considered to be invisible. The knowledge of "^o^le observed in a given spectral band to recover the bolometnc 

these mode visibilities is often requested for solar-like oscilla- ampiituae. 

tion analyses. They sensitively depend on the stellar limb dark- ^he fluctuations induced by radial modes identically affect 

the whole atmosphere, and can be considered as fluctuations 6T 

* Full Table 1 is only available in electronic form at the CDS via of the effective temperature T^g. For photomettic measurements, 

anonymous ftp to cdsarc . u- strasbg . fr ( 130 . 79 . 128 . 5 ) or via niade for example with Kepler, the oscillations generate fluctu- 

jhttp : //cdsweb . u- strasbg . fr/cgi-bin/qcat? J/A+A/5 3 1/A124 1 ations 6F}^ of the measured total flux F^. This flux depends on 
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Fig. 1. A solid line shows the spectral response &^ of Keple r 
as a function of the wavelength A (after IVan Cleve etal ] l2009h . 
whereas dotted and dashed lines show the spectral responses of 
CoRoT seismo- and exofield (Auvergneet al. 2009). Grey pro- 
files indicate the spectra (in arbitrary units) of black bodies with 
temperatures of 7000, 6000, and 5000 K, from top to bottom. 

the response of the instrument. By considering fluctuations small 
enough to be considered as linear perturbations, and by assum- 
ing that the observed star radiates as a black body of temperature 
Teff , the fluctuations 6Fy^ and 6T are linked through the relation: 



6Fk 



jTK(A)^(A,T,i,)6TdA 
jTK(A)B(A,T,n)dA 



(1) 



where B(A, T^g) is the Planck function, and Tk(A) is the transfer 
function of the instrument, and A the light wavelength (see M09 
for a detailed derivation). The transfer function reads 7k(A) - 
&Ki^)/Ev, with Sk the spectral response of the detector, and 
Ey - he I A the photon energy at the wavelength A. The constants 
h and c denote the Planck constant and t he light speed. The spec - 
tral response of Kepler is described in IVan Cleve et al.l (l2009lFl 
and plotted in Fig. [T] For comp arison purposes, the r esponses 
of the two CoRoT channels (see lAuvergne et al.ll2009l) are also 
plotted. 

The fluctuation 6Fb of the bolometric flux F\, is simply 
linked with 6T: 



Fb 

We then deduce 

SF\, SFk 

with the bolometric correction factor 
4 / TK(A)BiA,T,s)dA 



T,ff jTK{A)^{A,T,ff)dA 



(2) 



(3) 



(4) 



We numerically computed the bolometric corrections for ef- 
fective temperatures T^ff within the range 4000-7500 K, cover- 
ing the range of solar-like oscillating stars. Results are shown in 
Fig.|2] A second-order polynomial fairly fits the computed set of 
points. We then obtain the following law: 



CK-boKT'elf) = 1 + ai(Tef[ - To) + a2iT^s - To) 



(5) 



Fig. 2. Bolometric correction CK-boi as a function of the eff'ective 
temperature Teff . Crosses are numerical calculations of Eq.^] the 
line is a second-order polynomial fit of these points. 

with the coefficients To = 5934 K, oj = 1.349 x 10 K"', 
and fl2 = -3.12 x 10~' K"^. This approximation, also plotted 
in Fig. |2] gives consistent results with the numerical computa- 
tions within ~ 10 This error is probably negligible compared 
to the method error mainly induced by the departure of stellar 
spectra from black body emissions. 

Alternatively, we express the bolometric correction as a 
power law 



CK-hol(Tefi) - (Tefi/To)" 



(6) 



with a - 0.80. This approximation is consistent with the numeri- 
cal computations within ~ 10"^ over the considered temperature 
range. 

3. Visibilities of non-radial modes 

Computing the relative visibilities of non-radial modes com- 
pared to { = modes requires the knowledge of the stellar 
limb darkening. For a monochromatic observation at the wave- 
length A, the visibi Hty V[ of a rnode of degree { reads (e.g. 
iDziembowskill 1 9771: iOizon & Solankill2003l) 



Vc^ V^K2?TTj r Pe(M)W(fi)fidfi 
Jo 



(7) 



where P[ is the fth-order Legendre polynomial and W(jj) is a 
weighting function depending on fi, the distance to the limb 
(equal to at the limb and 1 at the disc centre), defined as the 
cosine of the angle MCO, where M is the considered point at the 
stellar surface, C the centre of the star, and O the observer. The 
function W{fi) is identified to the LD profile gAin) - h{p)l Lxi^), 
where Ixiii) is the specific intensity at the wavelength A. 

We also define the normalised -or relative- visibilities Vc - 
VtlV{), and the total visibility 



y2 

'tot 



phttp : //keplergo . arc . nasa . gov/CalibrationResponse . 
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We show (after lBallotll20T0[ Eq. A. 16) that 

/' 

Jo 



Vi = 2;r 



WQiff/dn. 



(8) 



(9) 



2 



J. Ballot et al.: Visibilities and bolometric corrections for stellar oscillation modes observed by Kepler (RN) 





Fig. 3. Total visibility V^^^ and relative visibilities for degrees i - 1,2, and 3 plotted as a function of the effective temperature 
Tetf , for different surface gravities log g (plotted with different line styles), and different metallicities [M/H] (plotted with different 
colours). 



Thus, we deduce the total visibility 
2 XV(//)Vd/z 



(10) 



For a broad band observation, the computation is not as sim- 
ple. In M09 an expression for t he mode visibilitie s is derived 
by following the approach of Ber thomieu & Provost! (il990l) . It is 
straightforward to reformulate their expression (M09, Eq. A. 18) 
in the same form as Eq. Qby considering W{fi) = Wk(m) with 



where 



(11) 



(12) 



Ga = J iJgA(fj)dfi and i/,, = ^J" g^(//) djuj . 

We computed LD profiles g,)(ju) as described in lBarban et al.l 
(l200l . For this, we used stellar atmo sphere models computed 
with the ATLAS9 co d^ (lKuruczl[T9"93i) in a modified version 



|http : //kurucz . harvard . edu | 



including the convective prescrip tion of Canuto et all (Il996l) . 
known as CGM (for details, see iHeiteret al.ll2002h . We then 
computed V^^^, as well as for £ - 1,2 and 3, for a large 
range of stellar fundamental parameters: for effective tempera- 
tures Teff € [4000,7500] K, surface gravities logg e [2.5,4.5], 
and three metallicitie^l [M/H] = -1.0,H-0.0, and H-1.0. Results 
are listed in Table [1] and are plotted in Fig. |3] As expected, it 
mainly depends on the temperature Tetf, while the sensitivity to 
log g is weak. However, strong changes in metallicity have visi- 
ble impacts on the results. 

In M09, a simplification has been done by noticing that 
the product H^Ga varies slowly with A. Moreover, V^^^ is not 
obtained through Eq. [TOl but computed as the truncated sum 
Zto As shown in Fig.|4]for logg = 4.5 and [M/H] = -hO.O, 
the results of this simplification give very close results to the 
more complete computations. 

Even by considering this simplification, the calculations of 
Vc and Vtot require one to know wavelength-dependent LD laws 
gAiM), of which computations are time-consuming. However, we 
find in the literature LD profiles which are already integrated 
over spectral bands of interest (for Kepler, see for example Sin^ 



3 We recall that [X/Y] means log(X/ F),„,-log(X/ Y)q, X and Y being 
chemical element abundances in number of atoms per unit volume. 
* This simplification was also made to obtain Eq.[T] see M09. 
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Table 1. Total visibility V^^^ and relative visibilities for de- 
grees £ - 1,2, and 3 for different fundamental stellar parameters 
Jeff, log ^, and [M/H]. 



Teif [K] 




[M/H] 


y'L 


vi 


Vi 




4600 


4.5 


-1.0 


3.17 


1.54 


0.58 


0.040 


4600 


4.5 


+0.0 


3.20 


1.55 


0.60 


0.044 


4600 


4.5 


+ 1.0 


3.23 


1.57 


0.61 


0.049 


5000 


4.5 


-1.0 


3.13 


1.53 


0.56 


0.034 


5000 


4.5 


+0.0 


3.16 


1.54 


0.58 


0.039 


5000 


4.5 


+ 1.0 


3.20 


1.56 


0.60 


0.044 


5400 


4.5 


-1.0 


3.09 


1.52 


0.54 


0.029 


5400 


4.5 


+0.0 


3.13 


1.53 


0.56 


0.033 


5400 


4.5 


+ 1.0 


3.17 


1.54 


0.58 


0.039 


5800 


4.5 


-1.0 


3.07 


1.51 


0.53 


0.027 


5800 


4.5 


+0.0 


3.10 


1.52 


0.54 


0.030 


5800 


4.5 


+ 1.0 


3.14 


1.53 


0.57 


0.035 


6200 


4.5 


-1.0 


3.06 


1.50 


0.52 


0.025 


6200 


4.5 


+0.0 


3.08 


1.51 


0.53 


0.027 


6200 


4.5 


+ 1.0 


3.12 


1.52 


0.55 


0.031 


6600 


4.5 


-1.0 


3.05 


1.50 


0.52 


0.024 


6600 


4.5 


+0.0 


3.06 


1.50 


0.52 


0.025 


6600 


4.5 


+ 1.0 


3.09 


1.52 


0.54 


0.028 


7000 


4.5 


-1.0 


3.05 


1.50 


0.52 


0.024 


7000 


4.5 


+0.0 


3.05 


1.50 


0.52 


0.024 


7000 


4.5 


+ 1.0 


3.07 


1.51 


0.53 


0.025 



Notes. The full table can be found at the CDS, or at 
^http://www.ast .obs-mip. fr/article998 .html| or at 
jhttp : / /www ■ j erome- ballot . f r/ , 

I2OIOI) . We denote g^ifi) = /K(y^)//K(l) the integrated LD profile 
where 

Ik(m)^ J TK(A)h(M)d^- (13) 

We propose to use Eq.|7]with W(fj.) - gKifA to perform approxi- 
mated computations of V[ and Vtot- Formally, this approximation 
is a correct simplification of the full computation for narrow- 
band filters only. Because the Kepler band is fairly broad, this 
approximation has to be verified first. We then computed the in- 
tegrated LD profiles gK for our different models, and used them 
to calculate V( and Viot- Results are shown in Fig.|4]for Vtot (Ve 
are not displayed, but their behaviours are the same as Vtot)- The 
obtained values for the different visibilities are satisfactory and 
reasonably close (better than 1 %) to the full computation. 

For comparison purp oses, we als o used the LD laws in the 
Kepler band published bv lSind (|2010() . We used its non-linear 4- 
coefficient LD lawfl to compute the visibilities. Results for Vtot 
are also shown in Fig. |4] (Here again, Ve show exactly the same 
behaviour). We then see noticeable differences with the results 
obtained with our own LD profiles, which are greater than the 
differe nce observed between our full and approximated compu- 
tations, has also made used of the ATLAS9 code for 
computing LD. The difference is mainly explained by different 
treatments of physical processes in the atmosphere models, es- 
pecially of the convection, as already pointed out in B arban et al, 
(l2003h . We recall that Sing (l2010l) used ATLAS9 model atmo- 
spheres computed for mixing-length theory with a mixing-length 
parameter a - 1.25 and overshooting included, while in the 
present work we used the CGM approach for the convection 

^ The improved 3-coefficient LD laws was also tested and yielded 
very close results. For a fair comparison, we also fitted our LD profiles 
with 4-coefficient laws and used them to compute visibilities: our results 
were almost unchanged. 



treatment without overshooting (see iHeiter et"an l2002[ for ar- 
gumentation). It shows the sensitivity of LD computations to the 
models and the consequences for mode-visibility calculations. 




3.0 r 



4000 4500 5000 5500 6000 6500 7000 7500 

T-effM 

Fig. 4. Total visibility V^^^ plotted as a function of the effective 
temperature Tefj, for logg - 4.5 and [M/H] - +0.0. The solid 
line shows the complete computation, the dotted-dashed line cor- 
responds to the simplification of M09, the dashed line shows the 
approximated results obtained by using integrated LD pro- 
files. The red dashed line is obtained by using the LD laws pub- 
Hshed by Sing (2010) . 



4. Discussion 

The notations used in this note differ from those in M09. The 
term CK-boi of the present work is equal to 4-/Rg in M09, and 
Vtot to Rose / Re ■ It is worth noticing that V^^^ is also equal to the 
c factor of iKieldsen et all (2008). Our results could appear to 
be more dependent on the metallicity than the computations of 
M09. This is simply an impression owing to the scales and the 
considered variables. For example, around Teff = 6000 K, V^^^ ~ 
3.1 and the spread of the different curves is around 0.06, then 
the corresponding dispersion for Rose is around 0.07, which is 
consistent with the point spread observed in Fig. 7 of M09. 

The coefficients introduced here appear quite naturally with 
the classical analysis techniques: 

1 . Using the factor CK-boi is straightforward: by multiplying the 
amplitude (for example in ppm) of an f = mode observed 
in the Kepler band, one recovers its bolometric amplitude. 

2. The visibility can be used to fix a priori ratios of mode 
heights in power-spectrum fitting, or to be a posteriori com- 
pared to the fitted ratios if they are free parameters. 

3. In some cases, specially for s tars with low signal-to-noise ra- 
tios (e.g. iMathur et al.ll2010h . the power density spectrum is 
convolved by a box car with a width equal to the large sepa- 
ration A. The smoothed spectrum is then used to recover the 
power maximum P„. Assuming that one mode of all degrees 
is present in the interval A and they all have a similar intrinsic 
amplitude (hypothesis of energy equipartition), we recover 
the maximum amplitude of radial modes in the Kepler band 
as Af=o,K = V-Pm/A/Vtot, and, then, the bolometric maximal 
radial-mode amplitude is Af=o,boi = CK-boi V^"™ /A/ Vtot- 
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At this point, it is important to stress the limits of the kind of 
approach we used. As already mentioned, our bolometric correc- 
tions ignore the departure of stellar surface emissions from black 
body spectra. As a consequence, the effects of photospheric ab- 
sorption lines and bands are neglected, as are the effects o f inter- 
stella r reddening. By using a classical reddening law (e.g lSchiidI 
the bolometric correction CK-boi is increased by less than 
2% for a reddening E{B -V) = 0. 14 -typical for the open cluster 
NGC681 1- and is negligible for nearby stars. 

Second, the LD laws are only models, which are also, in 
essence, approximations. It is nowadays possible to test mod- 
elled LD, at least linear approximations of LD, with various 
kinds of observations, such as interferometry (e .g. Burns et al.l 
[T997h . li ght curves o f eclipsing b inaries (e.g. iPopperl 1 19841) . 
planetary transits (e.g. 'Barge et al."2008"), or microlensing (e.g. 
[FouQue et al. 2010). Besides, Sing (2010) has shown discrepan- 
cies between its modelled LD and the LD deduced from plan- 
etary transits observed by CoRoT. Nevertheless, new 3-D at- 
mos phere modelling should provide better st ellar LD prediction 
(e.g.'Bi got et"al]|2006tlChiavassa et al.ll2009h . 

Last, even considering the LD profiles are perfectly known, 
some other assumptions could be not fulfilled. In our approach, 
the physics of modes inside the atmosphere is simplified: partic- 
ularly, non-adiabatic effects are neglected, and the stellar photo- 
sphere is assumed to be thin enough to render the variations of 
mode amplitudes with the altitude in the atmosphere negligible. 
For the Sun, whose surface is resolved, we have a direct mea- 
surement of the LD at different wavelengths. It is then possible 
to compute visibilities Ve from the real LD and to compare them 
to the ratios of observed mode amplitudes: Salabert et al.l (1201 ih 
have shown differences between the computed visibilities V( and 
the observations of the helioseismic VIRGO (Variability of Solar 
IRradiance and Gravity Oscillation) instrument. The lack of the 
previously mentioned physical effects in the modelled visibilities 
can be made responsible to explain the discrepancies. However, 
it is also possible that the intrinsic amplitudes of modes are dif- 
ferent. Other stars have show n unexpected high app arent am- 
plitudes of f - 3 modes (e.g. iDeheuvels etalj|2010 '). It is still 
unclear whether this is a consequence of incorrect computations 
of mode visibilities, or a signature of the over-excitation of some 
modes. These two examples show us how important it is to im- 
prove in the near future the computation of mode visibilities, by 
considering more realistic physics of acoustic modes in the stel- 
lar atmospheres, to be able to recover intrinsic mode amplitudes. 
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